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We present a detailed analysis together with exact numerical calculations on one-loop contribu- 
tions to the branching ratio of the radiative decay of fx and r , namely ^ — > e 7, r — > e 7, and r — > /.17 
from supersymmetry without R parity, focusing on contributions involving bilinear couplings. A 
numerical study is performed to obtain explicit bounds on the parameters under the present exper- 
imental limit. We present, and use in the calculation, formulas for exact mass eigenstate effective 
couplings. In this sense, we present an exact analysis free from approximation for the first time. 
After comparing our results against the closest early analysis, we discovered a major difference in 
resulted constraints on some [1* Bj combinations. Constraints from neutrino masses on the param- 
eters were considered. Our result indicates that the branching ratio measurement on /i — > e 7 down 
to 10~ 13 — 10 -14 and beyond, as targeted by the MEG experiment, has a chance of observing decay 
from the R-parity violating scenario. 
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I. INTRODUCTION 

Recent neutrino experiments have demonstrated that neutrinos change flavor as they travel from source to detector, 
a phenomenon consistent with the hypothesis of neutrino oscillation. All that contributes evidence for neutrino masses 
and lcpton-flavor violation (LFV) and provides the first definite experimental clue for physics beyond the standard 
model (SM). Many extensions of the SM predict a certain amount of LFV in relation to neutrino mass generation or 
otherwise. Important criterion for a viable model is giving acceptable neutrino mass spectrum while staying within 
experimental limits of LFV. Apart from the soft terms within the minimal supersymmetric standard model (MSSM), 
both LFV and neutrino masses are indeed forbidden by ad hoc discrete symmetry — the R parity. Note that soft terms 
by themselves still conserve total lepton number, and hence do not generate neutrino masses. In the supersymmetric 
standard model without R parity imposed, there is however an important source of LFV and neutrino masses. A 
major part of this comes simply from the R-parity violating (RPV) terms in the superpotential, though RPV soft 
(supersymmetry breaking) terms are also of interest. The latter too often escapes notice. 

The best evidence of supersymmetry (SUSY) would obviously be the discovery of SUSY particles in the collider 
machines. However, processes such as the leptonic radiative decays can serve as alternative ways to test SUSY, 
complementary to the direct SUSY particle searches. Although these processes have not yet been seen so far in 
present experiments, there are very stringent upper bounds on their possible rates implying important constraints on 
the new physics contributions. The present experimental upper bounds of branching ratio for r— ► [1] , r — *■ e7 
[2] , and /i — ► e 7 [3] are 

Bt{t -> /ry) < 6.8 x 1CT 8 , 
Br[r -> ej) < 1.1 x 1CT 7 , 
Br{p ej) < 1.2 x KT 11 . 

The muon radiative decay reaction /x — > e 7 has been the focus of most attention due to the experimental bound being 
much stronger. This bound will probably be improved in the future. The MEG experiment Ref. [4], which searches 
for [i — ► e7 decays down to 1CP 13 — 1(U 14 branching ratio is now in its final stage of preparation. The r decays may 
also be better probed in future facilities. 

The recent studies on radiative decays from other models such as little Higgs with T parity [5, 6] and the SUSY 
grand unified theories (GUT) [7, 8] also give some interesting results on the lepton-flavor violation processes. For 
example, in the little Higgs models with T parity, the presence of new flavor violating interactions and mirror leptons 
containing masses of order ITeV can enhance these processes to the level of the present experimental limit. In the 
SUSY GUT model, Ref. [7] discusses the complementarity between lcpton-flavor violation and Large Hadron Collider 
(LHC) experiments in probing the SUSY GUT. They found that the LFV experiments have strong capabilities to 
detect SUSY induced LFV, in some cases even outreaching the LHC. In Ref. [8], the authors study the correlation 
between J7 e 3 and the Br(/.i — > ej) in the context of a SUSY SO(10) framework. They find that taking running effects 
into account leads to a constant enhancement of the value of U e s at the high scale, bringing (i — > ej into the realm 
of MEG for SUSY parameter space regions which were previously excluded without inclusion of such running. 
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The study we present in this paper analyzes branching ratios that can be generated for all processes in the context 
of the generic supersymmetric standard model (GSSM), i.e. SUSY without R parity [9]. If one simply takes the 
minimal supersymmetric field spectrum of the SM and imposes nothing more than gauge symmetries while admitting 
soft SUSY breaking, the generic supersymmetric standard model would be obtained. Thus, the GSSM is the complete 
theory of SUSY without R parity, where all kinds of RPV terms are admitted without bias. Assuming SUSY, it 
is at least conceptually, the simplest model to accommodate neutrino mixing and oscillations. We work within the 
framework of single- VEV parametrization (SVP)[9-11], which is an optimal choice of flavor basis that helps guarantee 
a consistent and unambiguous treatment of all kinds of admissible RPV couplings and to maintain a simple structure 
for RPV effects on tree-level mass matrices for all states including scalars and fermions. 

The experimental bound on the branching ratio of these leptonic radiative decays is used to constrain the model 
parameter space, particularly the RPV part. Under constraints by the present experimental upper bound of branching 
ratios for fj, — > e 7 , t — > e 7 and r — > [i 7, we obtain the allowed region of the RPV parameter spaces. We give complete 
one-loop formulas for the type of contributions to the branching ratio of three leptonic radiative decays. We present 
numerical analysis of these contributions from all possible combinations of RPV parameters. Besides the more familiar 
Mfe ^kij 1 1 there are a list of combinations of type B"^ X^ij , [J,* Hj , and B* fij . A similar analysis on the /1 — + e 7 process 
has been reported in Ref.[13] in 2001 2 . The present work differs from Rcf.[13] in a few important ways. The present 
work is based on using new formulas of exact mass eigenstate couplings to calculate the one-loop diagrams, while in 
Rcf. [13] the authors only used electroweak states (7^'s) as an approximation for physical particles of external legs to 
the loop (the charge leptons). The latter amounts to neglecting the Higgsino and wino components of the decaying 
and product charged leptons. Therefore, the current analysis is an improvement or completion of the work reported in 
Rcf. [13]. In particular, we find that the constraints one can obtain on some of the ^1* Bj type parameter combinations 
have very substantial improvement. This is indeed the first exact calculation of processes within the model at the 
one-loop level presented, free from any approximation of the type. In addition, our results on the other two, r decays, 
processes have not been available in previous literature. We also increase the value of the fx parameter (corresponding 
to the MSSM /j, parameter) used from 100 to 135 GeV to accommodate the updated lightest chargino (xf ) lower mass 
limit of approximately 104 GeV [19]. In the sense explained above, the paper is somewhat of a sequel to Ref.[13], 
where we draw comparison when relevant. However, it can also be read just on its own. Readers who want to do so 
may simply neglect statements matching analysis and results here with that of the latter. One catch though is that 
we focus our discussion and result presentations on interesting results we get beyond that of Ref. [13], only briefly 
summarizing features and results that are essentially well explored in the latter reference. 

As experimental evidence for neutrino masses has become quite well established, we also include in our analysis 
a brief comparison of results from radiative decays with neutrino mass bounds. Recombinations of RPV couplings 
typically contribute both to the decays and neutrino masses, but with different dependence on the other model 



1 The interesting kind of combination of bilinear and trilinear RPV parameters contributing to flavor violations through scalar mass 
mixings [12] or a one-loop diagram [13, 14], and analogous one-loop dipole moment [15] were published a few years ago. More recently, 
similar contributions to radiative B decays were also published [16] . 

2 There have been numerous studies on similar processes from various versions or limited models of R-parity violation in the literature. 
Most of the model assumptions look ad hoc. We have no intention of reviewing all of that here. However, an early study on fi — » e 7 from 
softly broken R parity [17] should particularly be mentioned. Another particularly noteworthy paper on the topic is given by Ref. [18]. 
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parameter. Some RPV parameters, like the bilinear ones, can give rise to a neutrino mass term alone. However, we 
do not have solid evidence on the scale of the neutrino masses, only Am 2 . And there are so many relevant RPV 
combinations for both neutrino masses and radiative decay processes that it makes a comprehensive and systematic 
analysis unrealistic unless further assumptions are taken on the model structure. Some parameter combinations may 
have a more important role to play in a certain process while others may give dominant contributions to neutrino 
masses. We are interested in investigating and presenting generic results on model parameters. Hence, we adopt a 
naive strategy on the interpretation of neutrino mass bounds as naive upper bounds on the involved parameter(s) in 
order not to give a neutrino mass term contribution beyond the sub-eV scale. This may be a bit on the conservative 
side but is considered a reasonable strategy to be adopted. Note that this rough neutrino mass scale is not expected to 
be pushed down, since the scale of Am 2 is known from the oscillation experiments. On the other hand, we have only 
upper bounds for the radiative decays which may be, and we believe should be, probed with a lot better precision in 
the future. Our radiative decay results will hence be useful references for the future, even if they are no better than 
the naive neutrino mass bounds. In the case of fi — > e 7, for example, if the MEG experiment pushes the bound on 
the branching ratio down to 10~ 13 — 10~ 14 , our numerical results show that it can give a stronger constraint on the 
jU* B-2 parameter combination. To put it more interestingly, the neutrino mass bounds do not rule out the possibility 
of seeing a /1 -> e 7 at MEG coming from fx* B2 of the RPV supersymmetric model. In the case of (^2X212, the current 
bound is actually already competitive. 

This paper is organized as follows. In the next section, we briefly summarize the main features of GSSM and 
also set our notation. In Sec. Ill, we give the general exact formulas in the basis of mass eigenstates without 
any approximation. The focus is on the Cj — > £~ 7 amplitude from one-loop diagrams without colored intermediate 
states. The two sections are included here to make the paper self-contained, and set the notation to be used for the 
discussions that follow. Note that Sec. Ill does include important results, expressions for effective coupling among 
mass eigenstates involved, that have not been published before. Our numerical results will be presented in Sec. IV. 
We also compare the results obtained by using the exact mass eigenstate couplings versus the previous calculations 
and discuss the sources of the difference between them. In addition, we illustrate the effects of varying the input 
parameters on the bounds. Finally, Sec. V will be devoted to the conclusions. 



II. THE GENERIC SUPERSYMMETRIC STANDARD MODEL 



We briefly describe the model here. Details of the formulation adopted are elaborated on in Ref. [9]. The most 
general renormalizablc superpotential with the spectrum of minimal superficlds can be written as 



W=e ab 



H a H°L a + h% k Q1H u U£ + A'., ,,,. I ■'.',(.) ,P',, + - \ a f3 k LaLpE% 



2 XljkUiDjD% , 



(1) 



where (a, b) are SU(2) indices, and (i, j, k) are the usual family (flavor) indices (going from 1 to 3). The (a,/3) indices 
are extended flavor indices going from to 3. Note that A is antisymmetric in the first two indices, as required by 
the SU(2) product rules, shown explicitly here with e 12 = —s 2 i — 1. Similarly, A" is antisymmetric in the last two 
indices from SU(3) C , though color contents are not shown here. Besides the superpotential, the Lagrangian contains 
the gauge interaction part, including kinetic terms of the matter superfields and a soft SUSY breaking part. 

We take a definite flavor basis to write the model Lagrangian. Such choice of parametrization is not unique. In 
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the current case of the GSSM, the scalar parts of the colorless electroweak doublet supcrficlds could bear vacuum 
expectation values (VEVs). We use a parametrization called the SVP advocated by our group since Ref. [10]. A 
flavor basis with only one among the L a 's, designated as Lq, bearing a nonzero VEV is adopted. That is to say, the 
direction of the VEV, or the Higgs field Hd, is singled out in the four-dimensional vector space spanned by the L a 's. 
Explicitly, under the SVP, flavor bases are chosen such that (1) (Li) = 0, which implies Lq = Hd; (2) y| fc (= Xojk = 
- x J0k) = ^diag{m 1 ,m 2 ,m 3 }; (3) yf k (= X 0jk ) = & diag{m d , m s , m b }; and (4) y^ k 
where v = \* r 2(L ) and v u = \/2(H u ). 

The soft SUSY breaking part of the Lagrangian can be written as follows [9, 12] : 



V2 



*okm diag{ m u ,m c ,m t }, 



ae Tja f b fpC 



tab 



c 

3~k 



\A\ jk LfL)E° k 



h.c. 



-Q}m 2 Q Q + U^nfi U + D^ihl D + i^m 2 L L + &fhl E ■ 



+ h.c. 

t \H„ 1 2 



Mi ~ ~ 
~2 



fww 



^3 ~~ 

T99 



(2) 



where we have used Hd in the place of the equivalent Lo among the trilincar A terms. Note that iJmfL. unlike the 
other soft mass terms, is given by a 4 x 4 matrix. Compared to the MSSM case, mf 0CI corresponds to while f?if ofc 's 



give new mass mixings. 



III. LEPTONIC RADIATIVE DECAYS 



Within the GSSM, the three SM charged leptons are the light mass eigenstates out of a 5 x 5 charged fermions 
mass matrix, which also includes the charginos. We use the common notation n = 1 to 5, with the former states 
given by if = xt+2> i = 1 to 3. The states have characters different from the fermionic components l^'s and i+'s of 
the Li and Ef superfields, respectively, as a result of the generally nonzero fii RPV mixings between the charged 
leptons and charginos of the R-parity conserving (MSSM) limit. The smallness of the fa values as indicated by the 
resulted neutrino mass value [10, 20] was the basis for most of the approximations on related subject matters in the 
literature, essentially neglecting the difference between ^j's and k's. Reference [13] is not totally free from the kind 
of approximation, though it focuses on the //* A^fe RPV contributions to fj, — ► e 7. We will see that in this kind of 
parameter combination the approximation in Ref. [13], which neglects Higgsino and wino components of the decaying 
and product charged leptons, is perfectly fine. The study is the first of its kind, catching a major role of the /Xj's, as 
well as the soft bilinear RPV parameters -Bi's, in the LFV process in conjunction with the A-type couplings. 

As advertised, our analysis here goes beyond that. Let us start by looking into the full mass eigenstate couplings 
of the truly physical charged leptons. 



3 See, however, studies on the /i* ., RPV contributions for similar processes in the quark sector[15]. In fact, the general relevancy of 
the kind of parameter combinations to flavor diagonal and off-diagonal dipolc moments for fermions was first pointed out in Ref. [12]. 
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A. Charged scalar vertices 



A charged lepton if (or a generic Xn) couples to a charged scalar and a neutral fermion. From carefully expanding 
the Lagrangian, we have the vertices 



.92 



AC, 



■AC* 



*(Xn)0m +h-C. 



(3) 



where ^(1 =F 7s) ar c the L- and i?-handed projections and 
AC?™™ = U* n Xi n T> l 2m + U* n X^ k+4:)n V 



(k+2)r 



9i 



T> 1 



k+4)n ^(k+5)m 



1 

71 



[/^[tan^X^ + X/JPf^ 



1 

"71 

X 



U U+2)n [ tan ^-- X ln + X 2n] V 



(J+2) 



Cj'+2)n X 4n ^(j+5)m 



52 



WT -X", 



(j+2)n vl (fc+4))i ^(;i+5)m > 



(4) 



AC 



Ws ^4*i ^2m — ^Tn -X(fc+4)n ^ 



fc+4)n ^(fc+2)m 



1 

-V2 tan^V f * +2lfi X lra X>f 



77= [- tan6» vr X ln + X 2n ] T>[ r 



(j+2)n X (]+4)n u 2m 



-=*-v, 



92 



(j+2)n 



V' 



^khj 



v, 



0+2)r 



(j+2)n ^(k+4)n ^{h+2)m ' 



(5) 



ff2 W ' 32 

with n runs from 1 to 5, n from 1 to 7, and m from 1 to 8. V and 17 arc unitary matrices used to diagonalizc 
the charged fermion mass matrix, and X is the diagonalizing matrix of the neutral fermion mass matrix. T> 1 is the 
diagonalizing matrix of the charged scalar mass matrix [9]. We quote the corresponding terms A/^" m and Af t L nm from 
an earlier formula in Ref.[13] for comparison: 



ML 



i 

71 



[tan6WA7„ + X£ n ] V l (l+2)m 



= -V2 tanOwXinV 1 



(i+5)m 



92 



92 

XnD 1 



V* T) 1 

' ^4n U (i+S)m 



kih 

92 



X yk+4)n ^(/i+5)r 



(i+2)m 



Vet -^t- rr\l ^khi -i?- 
-\k+4)n L>2m X.k+4)n L>( 

92 <?2 



\h+2)m 



(0) 



One-loop diagrams formed by the pair of coupling vertices give rise to a class of contributions to the radiative decays 
we call neutralinolike, which obviously does include the ones with the physical neutralinos among the fermions running 
inside the loop. The new contributions come from the first five terms of Eq.(4) and first four terms of Eq.(5), which 
are not seen in Ref.[13]. These new terms involve higssinos and winos on the external legs and are easy to understand. 
For example, in N* mn , the first term denotes the interaction with wino, neutral Higgsino, and charged Higgs. The 
second term is the interaction with wino, i-handed sleptons and the neutrino flavor states, while the third term 
describes the interaction with charged Higgsino, neutrino, and i?-handed sleptons. The last two terms denote the 
interactions with charged Higgsino, charged Higgs, and the bino and wino, respectively. The nonzero /x's do give the 
physical charged leptons some Higgsino or gaugino components. 
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B. Neutral scalar vertices 



Next, we come to the charginolike contributions. Here we have to pair up neutral scalar vertices: 



92 *(X«) 



*(x;)C +h.c, (7) 



where 



,„ - i Z>|J - t/ 2 * fi Vi„ [B 



24 . 1 _ 24 . 

+ g7 V( - j+2)n ^/f ^(7+2)m _ * ^(7+7)™] _ U U+2)n V (j+2)n ~j= \P£m ~ 1 ^m] 



v 2 \ u ln -L [p- + i ] - v*. [/ 2 „ _L K ; - i ] 



24 1 24 1 

~ V (j+2)n U 2n -j= P(j+2)m + 1 V U+7)m] ~ ~^ V <J+ 2 )* U U+ 2 ) n 



-| bhl V* j+2 )n U( k +2)n — = \P(h+2)m + * ^(fc+rjm] ■ ( 9 ) 

with n and n run from 1 to 5 and m from 1 to 10. T> s is the diagonalizing matrix of the neutral scalar mass matrix 
[9]. Note that C^ mn is equal to C^ mn by definition. These are to replace Cf nm and C"* ln of Rcf.[13]. 4 

Cfnm = — Vln —7= [D{i+2)m + * 2?(i+7)m] — ~ V(i+2)n —?=■ Km _ * Am] 
V2 <?2 



92 \J 1 



Cinm — —U'2n —j= K+2)m + * Al+7)m] ~ ~ t^(j + 2)n —7= [Am + * Am] 
92 92 

jy (fe + 2)ii [2?(fe + 2 )m + * £fh + 7)m] • (10) 

92 \/ 1 

Note that V s is actually real, though we are using V s * notation as if it is not. This is just a convention for tracing the 
LFV structure of the various contributions in our analytical discussions below. Here, in fact, the real difference between 
the T> s * and V s terms is given explicitly by the different signs between the corresponding scalar and pseudoscalar 
parts. Note that the y,. terms in the above expressions can be written together with the A terms using the \ a pk 
notation and the identification of y c . as A u. This common structure between Lq and the L^s is very useful in our 
discussions below. 



4 Recall that in the latter case, one distinguishes the "external" charged lepton, approximated by an from an "internal" charged 
fermion mass eigenstate. 
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C. The decay amplitude 



In applying the above interactions to the process (j{p) — ► i\ j{q), we can write the amplitude as 



T = e e*° Ui(p - q) 



1 - 7s , A n 1 + 7: 



Uj(p) , 



m^iaapq' ( — r ^ 2 9 

where e* = e*(g) is the polarization four- vector of the outgoing photon. The decay rate is then simply given by 



|^ 2 | 2 + |^| 



(11) 



(12) 



It is straightforward to calculate the contributions from one-loop diagrams with the effective interactions of Eqs.(3) 
and (7). The result for A L 2 (A\ = A l 2 \l^r) is given by 



2 8tt sin 2 6U Ml 



(i + 2)r. 



.AC 



F 2 



'A/ 2 o 



\-2)r, 



Af R * 



^3 



Ml \ 



1 



2)„ 



F 2 



'51 
M 2 



87rsin 2 ^ M s 2 



M 2 



m, 5 1 M? 



/?£ /-lit 



M- /M 

-^6 



M 2 



(13) 



where 



F 2 (x) = 



1 



(l-6x + 3a: 2 + 2x 3 -6x 2 In a) 



6(l-x) 4 

F 3(x) = jz 1 ,, (1 - .x 2 + 2z In a) , 



F 5 (x) 



1 

6 (1 - a;) 4 
1 



(2 + 3x-6x 2 +x 3 + 6x lnx) 



(-3 + 4a;-x 2 - 2 lnx) 



(l-x)3 

with summations over all physical fcrmion and scalar mass eigenstates as represented by the n and m indices assumed. 

The processes we discuss here violate lepton flavor while conserving the lepton number. Before going into the 
analysis, it is instructive to introduce the lepton- flavor numbers L ei L^, and L T to the superfields as one does to 
their corresponding components in the SM. The RPV parameters bear violations of the lepton-flavor numbers. It is 
obvious that in order to have a contribution to £~ — ► l\ 7, a term must reduce Lj and increase Li by exactly one 
unit while leaving the rest unchanged. For instance, p\ ^ 2 means increasing a L e and reducing a L M while leaving L T 
unchanged. This simple but useful rule serves as a countercheck of individual contributions discussed below. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 



In this section, we present the results we obtained by a careful numerical implementation of our jj, — > ej, t — > cy, 
and t — > /1 7 formula with explicit numerical diagonalization of all the mass matrices involved. We isolate various 
major contributions by singling out each of the corresponding RPV parameter combinations as only nonvanishing one 
at a time. The soft SUSY breaking contributions to R-parity conserving slepton mixings are set to zero ( i.e. mf , mj, 
and A E are set to be diagonal ). A basic set of typical values chosen for the input parameters are given in Table I. We 
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used this set of inputs unless otherwise specified in the results below. Summary of bounds on various combinations 
of two R-parity violating parameters is shown in Table II. Neutrino mass bounds are also put onto the plots. We 
perform numerical calculation of the neutrino mass contributions from various parameters involved based on formulas 
from Ref.[21]. The bounds are obtained naively by requiring that individual neutrino mass terms obtained be less 
than the sub-eV level. Note that each of these bounds has different dependence on the background model parameters 
and are also typically different from the leptonic decay contribution term plotted. Hence, the comparison has only a 
simple illustrative value. A more comprehensive cross analysis is, however, considered not appropriate without taking 
further assumption on some of the parameters involved. At the end, we also show some of the effects of varying these 
input parameters. 



A. The |/x*A| or |S*A| contributions 



We first look at the contributions with a (fJ.*X) or (fiX*) structure. The dominant terms do not involve the Higgsino 
or wino component of the decaying or the product charged leptons. They come from the third term of Cf imn and the 
only one A-coupling vertex of C^ mn , corresponding to the diagrams with the chirality flip on the internal fermion line. 
Indeed, the role of the /li's come in through the internal fermion line, as discussed in good detail in Ref. [13]. Take 
^2, for example, the dominating term comes from C^ l , mn C^l nn , where n' < h. We then have the real scalar part of 
the contribution, for example, proportional to 

n=lm=l \ s m / U' 2 

In the \i — > e7 case, we have n' = 3 and h = 4 and then use the relation U^ + 2)n ~ <%+2),n an d Vy/_|_ 2 w ~ fi(j'+2).n'- 
The expression (14) can be given by 

EEKn M- U {k+2)n F 6 ( Vi; n V( h+2)m -fL . (15) 

n=lm=l \ s m j ^ 2 

If the loop function Fq could be factored out from the double summation, we would have a Vj* n M.- U^+2)n summation 
over fcrmions and a T^im^(h+2)m summ atk):ri over (real) scalars. Taking h = 2 in the above expression (15), we have 
the two dominating chargino contributions, the n = 1 and 2 parts, given approximately by 

\ ( M 2 _\ 

Kn^T 1 ^^ , (16) 



92 \ M £ 



s„ 



where R R is a 2 x 2 matrix with order 1 matrix elements (see Ref. [9] for details.) The expected combination pt* k \kn 
comes up, with k = 1 and 3 admissible. The same situation goes for the C^ mn C\ mn part of with the combination 
A£ 12 ( k = 2 and 3 admissible here ) instead. 

Terms involving the Higgsino and wino components of the decaying and product charged leptons are expected to 
be proportional to U afi or V an for a = 1 or 2, with [9] 

U a ( i+2 ) CC fa , 

Va(i+2) « Mi m % ■ ( 17 ) 



10 



In the case of the fj, — > e7, if we focus on the /Z3 and A 32 i parameter combination, the new contributions of A/" 3 ^, 
Af^ m ^, C^ n , and Cf^* do not contribute to the Br(/i — > ey), since they are only proportional to the \i\ or /12, not //3. 
Likewise, the /X3 A* 12 also does not have a new contribution to Br(fi — > cy). For the same reason, the contributions of 
//1 A* 23 , /Ltj A 132 , [ii AJ^, /X2 A231, and all B* A and B A* type combinations to the branching ratio of leptonic radiative 
decays are still essentially the same as the results obtained by neglecting the wino and higgino components of the 
decaying and product charged leptons. 

There are, however, new contributions that come from the charginolikc loop diagrams with the chirality flip on the 
external fermion line. They are the product of the fourth term and the A-coupling term of C? mn . There are two types 
of diagrams. The first type comes from C?/ mn C^ ln , where n' < fi, with a A coupling in the C r f mn and a Yukawa 
coupling in the C£, mn . It is given by 

ZJ U^n' U U+2)n V(f+2)n V* {k+2)n F 5 { j^f- J 2?y*+2) ro V {h+2)m " ~ ~ > ( 18 ) 



n—l m—1 



requiring further h = j' and k = f . Note that all off-diagonal matrix elements of the form V^+2)n are very small, 
those RPV ones (n — l or 2), in particular, contain a j/ e( , suppression [9]. A similar situation goes with the scalar 
sum 2^* +2 ) m 'D(h+2)m = fij'h by unitarity. m e stands for the mass of the decaying lepton. There is a factor of y e i 
suppression in the expression (18). The electroweak state Feynman diagram of this expression is in Fig. 1. Note the 
Higgsino component of an external line illustrated. 

In the case of r — + /17, we have h' = 4 and n = 5. After we take j' = h = k = 2 in the expression (18) and use the 
relation U(j + 2)n ~ <%+2),« an d U 2 \ i+2 ) a /*,■ , the expression (18) can be given approximately by 

Ml \ 

X, 

g 2 g 2 ~°\W 



* A232 Ve 2 T 

W -M 2 — — ^5 777" ■ (19) 



The expected u* A 232 comes out with a factor of y, 9 (muon Yukawa coupling) suppression. The same situation goes 
for the case of the fi — > ej and t —> ej, corresponding to the combination of //* A 121 and fi* A 131 , respectively, with a 
J4j suppression. 

Likewise, the second type comes from C^, mn C^" mn , where n! < n, with a A coupling in the C?, mn and a Yukawa 
coupling in the C^ n . It is given by 

U 2 n U* {j > +2 )n> ^{3+2)71 ^(fc+2)« F 5 \ Jft ) T> {j+2)m' D (h+2)m ^ ~^ ' ^ 



n—l m—1 



The electroweak state Feynman diagram for this case is in Fig. 2. In the case of a — > ej, we have h' = 3 and n = 4. 
After taking j = h = k = 2 and using the relation mentioned above, the expected /x 2 A* 12 would come out, with a 
factor of muon Yukawa. The same situation goes for the case of the r — > ej, and r — > ^7, corresponding to the 
combination of /i 3 A* 13 and /i 3 A* 23 , respectively, with a tau Yukawa. Compared with the first type, the second one 
does have larger contributions, since the former have the stronger Yukawa suppression. This result is confirmed by our 
exact numerical calculation and could be understood easily through Fig. 1 and Fig. 2, where the relevant illustrative 
electroweak state one-loop diagrams are given. Because these two types of contributions involve leptonic Yukawa 
couplings, they have Tjjrg dependence. However, they are not the dominant contributions, so the fi* A or fiX* type 
contributions are still insensitive to the tan/3. 
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We plot contours of the resulting branching ratio as a function of (real) fi 2 and A 2 i 2 , M3 and A 323 , and fi 3 and A 3 i 3 
in Fig. 3, Fig. 4, and Fig. 5, respectively. The present experimental limit is also shown and the allowed region at 90% 
C.L. is shaded. The 90% C.L. upper limit on Afc 2 i| or \/j,k A£ 12 | (normalized by |/x | = 135 GeV) is given by 

Im 3 ^32i I l^t ^i2i| 1^3 A* 12 | 1^2 A* 12 | <• 2 i x 10~ 7 (21) 

Im<>| ' |mo| ' ImoI ' |mo| 

Likewise, the 90% C.L. upper limit on |/i£Afc 32 | or |MfcA£ 23 | (normalized by | /x 1 = 135 GeV) is given by 

1/4 ^232! [Mi -^132 1 |Ma A* 23 | Imi A* 23 | <- 7 x 10~ 4 (22) 
lMo| ' I Mo I ' lMo| ' I Mo I 

The 90% C.L. upper limit on \n* k Afc 31 | or A£ 13 | (normalized by |/i | = 135 GeV) is given by 

1/4 -W| 1/4 ^131 1 I Ma A 313 | I Ma A* 13 | <- g 5 x 10~ 4 (23) 

Imo| ' Imo| ' |mo| ' |mo| 

For the BX* structure, one of the contributions is from the chirality flip inside the loop. In A L 2l this term comes 
from the A coupling term of C^, mn and the fifth term of C£ mn , where n' < n. We then have the expression 

EE V (j' + 2)fi' U (j+2)n V *U+2)n M x - n U(k+2)n F 6 [ 
m n=X 



■ [V* m + i^J [V( h+2)m + i% +7)m ] V f , (24 ) 

</2 i/2 

where the notation means the unphysical Goldstone mode is omitted. In the case of fj, — > ej, we have nf = 3 
and n = 4. Using VTji+2)n' ~ <5(j'+ 2 )™' ano - ^b'+ 2 )« ~ &(j+2)n an d taking k = 2, we then have the expression 



5 /M 2 _ 



£ £ K„ M ^4„ F 6 ( ] [2? 2 s m + iVU [V { s h+2)m + iV( h+7)m ] ^ . (25) 

The fermionic sum suggests a major contribution from n = 4, i.e. the muon itself with the m M dependence. However, 
this contribution then has two factors of "muon Yukawa" (y. 2 ) suppression. The scalar sum gives a contribution 
proportional to B*, tan/3. The expected combination B* h A^a comes up with an explicit tan [3 dependence. The other 
contribution having similar strength to the one in expression (24) discussed above is from the chirality flip on the 
external muon line. In An, this term comes from the A coupling term of C-, and the fifth term of Ck mn , where 

z; f o n' ran nmn' 

n 1 < n. We then have the contribution proportional to 

5 5 ( M x- \ » -Xhkj> -V>j 

E E V U+2)fi V *(]'+2)n> Uy+^n U (k+2)n F 5 \ —jg- ^>( h+ 2) m V 2m ~ ^~ — ~ • ( 26 ) 



m , 

n—l m— 1 



M? <? 2 <? 2 



In the case of \i — ► e7, we have n' = 3 and n = 4. Using V(j + 2)fj ~ <5(j+2).n and Vqv + 2W ~ \j'+2)n' and taking 
k = 2, we have the expression 

5 5 / Ml 



™» E E ut. tf*. j%( ^ j m + 2 )m vr m , (27) 



where we have ft = 3 or 1. The scalar sum gives the dominating contribution proportional to B* h . The expected 
combination B* h \ nTL comes up. The 90% C.L. upper limit on \B* k Afc 21 | or \Bk A£ 12 | (normalized by |m | = 135 GeV) is 
given by 

1^3* A321I \B* A 121 | \B 3 A 312 | \B 2 X* 212 \ ^ 1qv1n -4 / 08 n 
iMol 2 ' iMol 2 ' iMol 2 ' iMol 2 <L3X10 ■ ^ 
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Likewise, the 90% C.L. upper limit on \B£ Afe 32 | or \B k A£ 23 | (normalized by \fi \ = 135 GeV) is given by 

\ B 2 ^232 1 \B* A132J \Ba A,3 23 l \Bi A* 23 | 1n _3 /9Q N 

iMol 2 ' iMol 2 ' U 2 ' iMol 2 <L4Xl ° ■ W 

The 90% C.L. upper limit on \B* k A fc31 | or \B k A£ 13 | (normalized by |/x | = 135 GeV) is given by 

t^T' t^F' "ki 5- ' ~w <L9x10 ■ (30) 

B. The |-B*/x| contributions 

Next, we will discuss the B* fj, and B fx* type combinations. The dominant terms come from two types of diagrams 
with the chirality flip inside the loop. One of them comes from the third term of C£ mn and the fifth term of C^ mn . 
In A 2 , this term comes from C?, mn C-nmni where ft! < n. We then have the real scalar part of the contribution 
proportional to 



5 5 / m 2 



E E E7 O v +2)n' V (j+2)n V ln Af -U*J+2) n Fq\ jjf- J Py, +2)m - 1 ■ (31) 

n=lm=l V Sm / y2 

The relevant electroweak state Feynman diagram is then in Fig. 6. With only l i on the external legs, it obviously 
involves no Higgsino or wino component there. The dominating part with the charginos (n = 1 and 2) gives a 
fij dependence through U*(j +2 \ n . The scalar sum gives the contribution proportional to B*, tan/3. The expected 
combination B*, \ij comes up with an explicit tan/3 dependence. Similarly, in A 2 , this term comes from C^, mn C£ mn , 
again n' < h. It is given by 



5 5 / M 2 



E E V (0'+2)n' U (j+2)n U(j> +2 )n M x ;^ln F e[ TpT ) V 2m V (j+l)m ~ ~ ■ ( 32 ) 
n=lm=l V s ™ / ^ 2 

The expected combination /i*, Bj comes up (see Fig. 7). 

The other type of contribution comes from the second term of C£ mn and the fourth term of C^ mn , which is a new 
contribution not considered in Rcf. [13]. It docs have a sizable contribution in some RPV parameter space regions. 
In A 2 , this term comes from C^, mn C^ mn . We then have the contribution proportional to 

JL / M 2 

'X-n 

w 



E E ~ U 2n' V (j+2)n Vln M x -U*2n F e[ TT^ ) V 2m V (j+2)m ~f- ■ ( 33 ) 



„=lm=l \"* s ™/ g ' 2 

We expect a Bj tan/3 from the scalar mixing part, and a from U" 2n i. In A 2 , this term comes from C^ mn C|^ 

It is given by 



5 5 / M 2_ 



E E ~V(f + 2)n> U ™ U 2n M- V* ln Fei ) ^, +2)m ^ . (34) 

n=lm=l \ Sm / ^ 2 

Again, we expect a -B*/ tan/3 from the scalar mixing part, and a /i fi _ 2 from the U 2n . The relevant electroweak state 
Feynman diagrams for the cases are given in Figs. 8 and 9, respectively. They clearly illustrate the role of the Higgsino 
component on an external line. 

Taking the /i — * ej for example, we have n 1 = 3 and n = A. The expression (31) would become 

5 5 / M - \ 

E E V ^ K-^n F 6 [ j^- Vi m VZ £l , (35) 



n— 1 m—1 



M2 I Am Am ^ 
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where we used the relation C^y/_|_ 2 )3 ~ an d V(j+2)4 ~ ^,2- As mentioned above, the expected -Bj* /-t 2 combination 
comes up. It clearly has a muon Yukawa suppression. Likewise, the expression (32) becomes 

5 5 / M - \ « 

n=lm=l \ s m y i/2 

It clearly has the fj,* B 2 combination with an electron Yukawa suppression, and is thus smaller than expression (35). 
It is confirmed by our numerical calculation. The result we mentioned above is also discussed in Ref. [13]. The 
expression (33) becomes 

E E -^23 V ln M - lT 2n fJ-^) Vi m VZ , (37) 

and we have the n\ B 2 combination with a muon Yukawa suppression. The expression (34) becomes 

55 ( m2 -\ v 

E E K„ K-Um F G [ -T^t ^im 9& — , (38) 

and we have the B\ fi 2 combination with an electron Yukawa suppression. We expect that this contribution is smaller 
than in expression (37), because it has a larger Yukawa suppression. 

In summary, using an approximate formula of mass eigenstate couplings, one can only obtain the expression (35) 
and (36). However, in this paper, we take all contributions into account including diagrams with higssinos and winos 
on the external legs. Therefore, additional contributions, such as the expression (37) and (38) in this case, would be 
obtained. If we consider the B\ fi 2 combination, the dominant term still comes from expression (35), since the new 
contribution, expression (38), has an electron Yukawa suppression, while expression (35) only has a muon Yukawa 
suppression. However, the situation would be totally changed in the combination of //J B 2 , since the expression (36) 
has an electron Yukawa suppression while the new contribution, expression (37), only has a muon Yukawa suppression. 
As a result, our exact formula made the allowed region of the (£?2, Mi) parameter space more stringent than that in 
Ref. [13]. We give contours of Br(/j, — > cy) in the real (i?2,Mi) plane in Fig. 10. The solid lines represent the results 
obtained by using the exact mass eigenstate couplings, while the dashed lines stand for the results obtained by using 
the approximate ones. Likewise, the contour of Br(r — > cy) is given in the real (-B3, /Ji) plane in Fig. 11. The present 
experimental limit is also shown and the allowed region at 90% C.L. is shaded. 

The 90% C.L. upper limit on |/_t* B\ or |/U-B*| type combinations (normalized by \fi \ 3 , |/i | = 135 GeV) is given by 



< 6.5x10"', ' , < 7.1 x 10 

ImoI 3 ImoI 3 



-7 



lB *^ 1 <1.4xl0- 4 , i^<1.5xl0- 4 



lMo| 3 ' ' l/zol 3 



< 1.1 x 10-*, L^<i. 2 xl0- 4 . (39) 



1AM I Mo 

For a direct contrast, we also give the incorrect upper bounds obtained by using the approximate formula: 

|B2M *' <1.4xHT 4 , ^i^< 3.5x10°, < 2.3 x 10- 3 . 



ImoI 3 ■ ' ImoI 3 ■ ' ImoI 3 
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We can see clearly that the ratio between the limit of and the correct one in Eq.(39) is approximately the ratio of 

^2-. Therefore, the difference of the experimental bounds between the exact formula (solid line) and the approximate 
formula (dashed line) in Fig. 10 is the result of different Yukawa suppressions, which is consistent with our analysis. 
If we take the expected improvement from the MEG experiment into account and assume Br(fi — > ej) < 10~ 14 , we 
get 

L^i< 1.2x10-8, ^n<i.3xl(T 8 . 
Notice that this constraint is even more stringent than the one from neutrino masses as you can see in Fig. 10. 

C. Parameter variations 

In this section, we illustrate the effects of varying the input SUSY parameters on the bounds, using A 121 1 and 
B 2 \ as examples. The results are summarized in Table III. In the table, we list the variation of the /i and the 
M 1 in parts i and ii, respectively Our numerical results show that the bound is most stringent for small |/i | in both 
A 121 | and \fi* B 2 \ cases. In addition, the increase of M x also weakens the bound. These results are reasonable since 
increasing /i and AI X = |M 2 essentially increases the chargino and neutralino masses. In the case of A I2I |, the 
dominant diagram involves mainly the 1°, while the |/x*-B 2 | case involves the mixing between 1° and 1°. Therefore, 
varying does not have much effect on the bounds while varying the corresponding entries in m£ changes the bounds 
significantly (see parts iii and iv). 

Finally, part v of Table III shows the tan/3 dependence of the results. From the table we can see that varying tan/3 
has only a little effect on |/x* X 121 \, but a rather significant effect on |/x* B 2 \. Although the subdominant contributions 
mentioned in Sec. IV B involve the Yukawa couplings, and thus have the co * „ dependence, the dominant contributions 
do not have the tan/3 dependence [13]. As a result, the lack of sensitivity to tan/3 in the former case is to be expected. 
In the latter case, the numerical result shows that the bound has a strong dependence on tan/3. There are two sources 
that result in this dependence. The first one is the dependence of the Yukawa coupling. The other is the explicit 
tan /3 dependence of the dominant terms mentioned in Sec. IV C. Figure 12 shows a contour plot of the experimental 
bound of B(fi — ► e 7) in the (real) plane of (B 2 , /Xi) for various values of tan/3. It not only shows the tan/3 dependence 
of the results, but also illustrates that the experimental bound gives a more stringent constraint in the large tan/3 
region. 

V. CONCLUSION 

In this paper, we have given explicit formulas and detailed discussions on the full one-loop contribution to the 
radiative decay of fi and r, namely, fi — ► e 7, r — > e 7, and r — > /X7 for the generic supersymmetric SM (without R 
parity) . We use the exact formula of the mass eigenstate couplings to calculate the branching ratios of these leptonic 
radiative decay processes and compare them with the results obtained in an earlier approximation in Ref. [13]. In 
some combinations of RPV parameters, the results obtained by using these two approaches are exactly the same such 
as fi± A* 23 , \i\ A132, /X2 A 213 , \i\ A 23 i, /I3 A321, /X3 A* 12 , and all B* X and B A* type combinations. The dominant terms of 
the other /x*A or /xA* combinations such as fi2 A 212 , fii A 121 , /i 3 A* 13 , fx* A 131 , /X3 A 323 , and fi 2 AJ 32 are still the same in 
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the two methods, but the subdominant terms come from a new contribution. In the B* \i or B /i* contributions, the 
dominant contributions of B* fj, 2 , B* /x 3 , and B* fi 3 are the same in the exact and approximate formulas. However, the 
dominant contributions of /i* B 2 , fx* B 3 , and /i* B 3 are totally new. The upper bound on these combinations obtained 
from the experimental limit are 

\^o\ 3 

< 1.5 x icr 4 , 

< 1.2 x icr 4 . 

As a result, our exact formulas impose more stringent constraints on the admissible region of parameter spaces for 
the GSSM, or SUSY without R parity. If we also consider the expected improvement from the MEG experiment and 
assume Br(fi — > ej) < 10~ 14 , we could get even more stringent constraints on the B* combination 

< 1.3 x icr 8 . 

The constraint is even more stringent than the naive constraint from neutrino masses imposed here, which indicates 
a very encouraging scenario for future probing of the leptonic radiative decay and may be as well as the r decays. A 
more involved analysis will have to be performed on the full model parameter space matching the radiative decays to 
neutrino mass generations. We want to highlight though that the scale for the actual neutrino masses is not expected 
to be reduced while probes on the leptonic radiative decays, and for that matter the other lepton number/flavor 
violating decays, can be improved. That makes the latter a promising ground to further explore models like SUSY 
without R parity with rich lepton number/flavor violating structures. 
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TABLE I: Basic input SUSY parameters for the numerical results presented. These values are adopted unless otherwise 
specified. 



Mi (GeV) 


M 2 (GeV) 


fio (GeV) tan/3 


100 


200 


135 40 


ml (10 4 GeV 2 ) 


mj (10 4 GeV 2 ) 


A e (GeV) 



diag{2, 1, 1, 1} diag{l,l,l} 100 
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FIG. 1: Diagram with chirality flip on the external line, which has larger Yukawa suppression. 
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FIG. 2: Diagram with chirality flip on the external line, which has smaller Yukawa suppression. 
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Br(^ --> e y ) 




Neutrino mass bounds 



1E-6 1E-5 1E-4 1E-3 0.01 0.1 1 10 

H 2 (GeV) 



FIG. 3: Contours of B(/i — > e'y) in the (real) plane of (/i2,A2i2). The 90% C.L. allowed region is shaded. The dash-dotted 
(red) lines are sub-eV neutrino mass bounds. 



TABLE II: Summary of bounds on various combinations of two R-parity violating parameters, normalized by j/i | = 135 GeV. 
The input parameters are as in Table I. 





M3 A321 1 


Ml ^12ll 


|M3 A312I 


or 


|M2 A 2 i2 1 


< 2.1 x 10" 


-7 




Imo ' 


MO 1 ' 


MO 1 ' 


MO 






M2 -^232 1 


Ml A132I 


|M3 ^323 


or 


Imi A123 1 


< 7.0 x 10" 


-4 




ImoI 


ImoI 


MO 1 ' 


MO 






M2 ^231 


Mi A131I 


|M3 A313 


or 


|M2 A|i 3 1 


< 8.5 x 10" 


-4 




ImoI 


Imo 


MO 1 ' 


MO 






B\ M2 

ImoI 3 










< 6.5 x 10" 


-7 




B 2 M*l 
Mol 3 










< 7.1 x 10" 


-7 




B l M3 

Mn 1 3 










< 1.4 x 10" 


-4 




S3 M*l 

MO I 3 










< 1.5 x 10" 


-4 




B 2 M3l 
MO 1 3 










< 1.1 x 10" 


-4 




S3 M2 1 

MO 1 










< 1.2 x 10" 


-4 




Mi M2 
Mol 2 










< 3.7 x 10" 


-5 




M^ M3 
M0 I 2 










< 4.7 x 10" 


-3 




M2 M3l 
1 2 
MOT 










< 3.6 x 10" 


-3 




S3 A321 | 


\B't A121I 


IS3 A3 12 | 


, or 


S 2 AJ 12 | 


< 1.3 x 10" 


-4 




ImoI 2 ' 


ImoI 2 


Imo I 2 


r> 

MO 






S| A232 1 


\B'l A132I 


1-^3 -^323 I 


, or 


Si A 193 


< 1.4 x 10" 


-3 




ImoP ' 


ImoI 2 


1 i2 

Imo r 


12 
MO 1 






S2 A231 1 


\B\ A131I 


1^3 A5 13 | 


, or 


l S 2 A213I 


< 1.9 x 10" 


-3 


ImoI 2 ' 


ImoI 2 


MO 1 


2 

MO 1 
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FIG. 4: Contours of B(t — > ^7) in the (real) plane of (/i3,A 32 3). The 90% C.L. allowed region is shaded. The dash-dotted 
(red) lines are sub-eV neutrino mass bounds. 
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0.01 



r> 1E-3 
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Neutrino mass bounds 



m (GeV) 



FIG. 5: Contours of B(t — > ey) in the (real) plane of (^3^313). The 90% C.L. allowed region is shaded. The dash-dotted 
(red) lines are sub-eV neutrino mass bounds. 
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FIG. 6: The charginolike loop diagram contributes to the leptonic radiative decay due to Bp fj,j combination. 
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FIG. 7: The charginolike loop diagram contributes to the leptonic radiative decay due to /x*, Bj combination. 
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FIG. 8: The charginolike loop diagram contributes to the leptonic radiative decay due to Bj /U*/_„ combination, where n! — 2 < j 
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FIG. 9: The charginolike loop diagram contributes to the leptonic radiative decay due to Bp combination, where j' < n — 2 
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Br(n -> e y) 




Exact 

Approxmation 
Neutrino mass bounds 



B 2 (GeV 2 ) 



FIG. 10: Contours of B(fi — > e 7) in the (real) plane of (B 2 ,fJ,i). The 90% C.L. allowed region is shaded. The solid lines represent 
the results obtained by using the exact mass eigenstate couplings, while the dashed lines stand for the results obtained by using 
the approximate ones. The dash-dotted (red) lines are neutrino mass bounds. Notice the MEG experiment targets probing the 
decay at 10~ 13 - 10 -14 . 



Br( t --> e y) 



=r 0.1 



0.01 





° ■■ 


Exact 






Approxmation 

Neutrino mass bounds 


: ^XX 

s 














1 I 1 1 1 





1 10 



B 3 (GeV 2 ) 



FIG. 11: Contours of B(r — -> e 7) in the (real) plane of (B 3 , ^x). The 90% C.L. allowed region is shaded. The solid lines represent 
the results obtained by using the exact mass eigenstate couplings, while the dashed lines stand for the results obtained by using 
the approximate ones. The dash-dotted (red) lines are neutrino mass bounds. 
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FIG. 12: Contours of the experimental bound of B(fi — » e*y) in the (real) plane of (B 2 , /Hi) with different values of tan/3, ranging 
from 10 to 50, among which tan/3=50 gives the most stringent constraint and the corresponding 90% C.L. allowed region is 
shaded. 
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TABLE III: Effects of parameter variations of interest, on the bounds of |m* A121 1 • (135 GeV) 1 and \fi* B 2 \ ■ (135 GeV) 3 . Note 
that the fixed mass scale of 135 GeV is used for normalization to extract numerical bounds. 

Parameter changes Normalized numerical bounds 

Mi >121 I Mi -B2I 
(135 GeV) (135 GcV) ? 

Original inputs of Table I < 2.1 x 1(T 7 < 7.1 x 1CT 7 

(i) ml = diag{20000, 500 2 , 500 2 , 500 2 } GeV 2 



ml =diag{500 2 ,500 2 ,500 2 } GeV 2 





Mo = 


- 500 GeV 


< 7.5 x 10" 


-6 


< 2.1 X 10" 


-4 




Mo = 


- 250 GeV 


< 2.7 x 10" 


-6 


< 3 X 10" 


■5 




Mo = 


— LOO Vjrc V 


*• ±.0 A 1U 


-6 


^ 71 y 111" 
1 . L A L\J 


■6 




Mo = 


- _1 S5 dpV 

— i-OO VjrC V 


< 1.3 x 10" 


-6 


< 7.3 x 10" 


■6 




Mo = 


- —250 GeV 




-6 


< 3.1 x 10" 


-5 




Mo = 


- —500 CleV 


/ so y in- 


-6 


< 2.1 x 10" 


-4 


(ii) Mi = 


1 71 .f _ cnn p„v 

2IVJ2 — JUU Vjrc V 


< ±.Z A 1U 


-6 


x J C y 1(1" 
<^ *±.0 A 1U 


-6 


(iii) 


m 2 = 


- ZUUUU X Qldg^l, 1, 1, Ij- Vjre V 


^ 9 q v 1 n~ 


-7 


Q 7 v in - 
\ y. ( x iu 


■7 




ml 


— Qiag^ZUUUU, 1UUU , 1UUU , lUUU j VjG V 


^ q n ^ 1 iv 

<^ O.U X 1U 


-6 


<. Z.l X 1U 


-5 


(iv) 


m% = 


- zUUUU x uiag|l, l, ij- Kjey 


<^ Z.l X 1U 


-7 


^ 7 o w i ri- 
< /.o X 1U 


■7 




mj 


= diag{1000 2 , 1000 2 , 1000 2 } GeV 2 


< 2.2 x 10" 


. 7 


< 8.7 x 10" 


. 7 


(v) 


ml = 


: diag{20000, 500 2 , 500 2 , 500 2 } GeV 2 












ml 


= diag{500 2 ,500 2 ,500 2 } GeV 2 












Mo 


135 GeV, tan/3=2 


< 8.7 x 10" 


- 7 


< 2.6 x 10" 


3 






tan/3=10 


< 1.1 x 10" 


-6 


< 1.1 x 10" 


-4 






tan/3=50 


< 1.3 x 10" 


-6 


< 4.5 x 10" 


■6 




/; = 


250 GeV, tan/3=2 


< 1.7 x 10" 


-6 


< 1.2 x 10" 


-2 






tan/3=10 


< 2.4 x 10" 


-6 


< 4.7 x 10" 


-4 






tan/3=50 


< 2.7 x 10" 


-6 


< 1.9 x 10" 


-5 




Mo = 


500 GeV, tan/3=2 


< 3.9 x 10" 


-6 


< 8.3 x 10" 


-2 






tan/3=10 


< 6.4 x 10" 


-6 


< 3.2 x 10" 


-3 






tan/3=50 


< 7.5 x 10" 


-6 


< 1.3 x 10" 


-4 



